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This Letter presents the first search for magnetic monopoles with the ATLAS detector at the CERN 
Large Hadron Collider using an integrated luminosity of 2.0 fb 1 of pp collisions recorded at a center- 
of-mass energy of As = 7 TeV. No event is found in the signal region, leading to an upper limit on 
the production cross section at 95% confidence level of 2 fb for Dirac magnetic monopoles with the 
minimum unit magnetic charge and with mass between 200 GeV and 1 500 GeV. No assumption about 
the production mechanism is made for this result, which is valid in the fiducial region given by pseudo- 
rapidity \r]\ < 1.37 and transverse kinetic energy 600-700 < E kin s'm8 < 1400 GeV. The minimum value 
of 700 GeV is for monopoles of mass 200 GeV, whereas the minimum value of 600 GeV is applica- 
ble for higher mass monopoles. Assuming the kinematic distributions from Drell-Yan pair production 
of spin-1/2 Dirac magnetic monopoles, the upper limit on the cross section at 95% confidence level 
varies from 145 fb to 16 fb for monopoles with mass between 200 GeV and 1200 GeV. 
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This Letter presents a search for magnetic monopoles with the ATLAS detector at the CERN 
Large Hadron Collider using an integrated luminosity of 2.0 fb -1 of pp collisions recorded at a center- 
of-mass energy of y/s = 7 TeV. No event is found in the signal region, leading to an upper limit 
on the production cross section at 95% confidence level of 2 fb for Dirac magnetic monopoles with 
the minimum unit magnetic charge and with mass between 200 GeV and 1500 GeV. No assumption 
about the production mechanism is made for this result, which is valid in the fiducial region given 
by pseudorapidity \t}\ < 1.37 and transverse kinetic energy 600-700 < E hla sm6 < 1400 GeV. The 
minimum value of 700 GeV is for monopoles of mass 200 GeV, whereas the minimum value of 
600 GeV is applicable for higher mass monopoles. Assuming the kinematic distributions from Drell- 
Yan pair production of spin-1/2 Dirac magnetic monopoles, the upper limit on the cross section at 
95% confidence level varies from 145 fb to 16 fb for monopoles with mass between 200 GeV and 
1200 GeV. 

PACS numbers: 14.80.Hv, 29.20.db, 29.40.Cs, 29.40.Vj 



Magnetic monopoles have long been the subject of ded- 
icated search efforts for three main reasons: their intro- 
duction into the theory of electromagnetism would re- 
store the symmetry between electricity and magnetism 
in Maxwell's equations; their existence would explain the 
quantization of electric charge pQ; and they appear in 
many grand unified theories However, to date no 
experimental evidence of a magnetically charged object 
exists. 

Recent searches for magnetic monopoles from astro- 
physical sources [SHS] are complemented by searches at 
colliders [10HI4] . This Letter describes a search for mag- 
netic monopoles in proton-proton collisions recorded at 
a center-of-mass energy of s/s = 7 TeV using the ATLAS 
detector at the CERN Large Hadron Collider (LHC). 

The Dirac quantization condition [T], given in Gaus- 
sian units, leads to a prediction for the minimum unit 
magnetic charge g: 
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where e is the unit electric charge and a e is the fine 
structure constant. With the introduction of a magnetic 
monopole, the duality of Maxwell's equations implies a 
magnetic coupling [15] 
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where /3 = v/c is the monopole velocity. For relativistic 
monopoles, a m is very large, precluding any perturbative 
calculation of monopole production processes. Therefore, 
the main result of this analysis is a fiducial cross-section 
limit for Dirac monopoles of magnetic charge g derived 
without assuming a particular production mechanism. 
A cross-section limit assuming the kinematic distribu- 
tions from Drell-Yan monopole pair production is also 
provided. 



Monopoles are highly ionizing particles, interacting 
with matter like an ion of electric charge 68. 5e, accord- 
ing to Eq. [T] The dominant energy loss mechanism for 
magnetic monopoles in the mass and energy range con- 
sidered herein is ionization |16H18j . The monopole ion- 
ization stopping power |18j , which is independent of the 
monopole speed /? to first order, results in the produc- 
tion of a large number of 5-rays. These energetic "knock- 
on" electrons emitted from the material carry away en- 
ergy from the monopole trajectory and further ionize the 
medium, fn the mass and energy regime of this study, 
the 5-rays have kinetic energies ranging from 1 MeV to 
a maximum of ~ 100 MeV. The secondary ionization by 
these (5-rays represents a significant fraction of the ion- 
ization energy loss of the magnetic monopole |18) . 

In the ATLAS detector [HO HO] , the monopole signa- 
ture can be easily distinguished using the transition radi- 
ation tracker (TRT) in the inner detector and the liquid 
argon (LAr) sampling electromagnetic (EM) calorimeter. 
The TRT is a straw-tube tracker that comprises a barrel 
< 1.0) with 4 mm diameter straws oriented parallel 
to the beam-line, and two endcaps (0.8 < \r/\ < 2.0) 
with straws orientated radially. Energy deposits in a 
TRT straw greater than 200 eV (low-threshold hits) are 
used for tracking, while those that exceed 6 keV (high- 
threshold hits) typically occur due to the transition radi- 
ation emitted by highly relativistic electrons when they 
penetrate the radiator layers between the straws. How- 
ever, the high-threshold hits can also indicate the pres- 
ence of a highly ionizing particle. A 2 T superconduct- 
ing solenoid magnet surrounds the inner detector. The 
LAr barrel EM calorimeter lies outside the solenoid in the 
|?7| < 1.5 region. It is divided into three shower-depth lay- 
ers and comprises accordion-shaped electrodes and lead 
absorbers. The cell granularity in the second layer is 
A?7 x A(f> — 0.025 x 0.025. The characteristic signature of 
magnetic monopoles in ATLAS is a large localized energy 
deposit in the LAr EM calorimeter (EM cluster) in con- 
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junction with a region of high ionization density in the 
TRT. A search for particles with large electric charge, 
which yield a similar signature, was performed previ- 
ously [21 and production cross-section limits for such 
particles were set |22j . 

The trajectory of an electrically neutral magnetic 
monopole in the inner detector is straight in the j — <f> plane 
and curved in r—z. The behavior of magnetic monopoles 
in the ATLAS detector is described by a Geant4 [23] 
simulation |24j , which includes the equations of motion, 
the ionization, the 5-ray production and a modified Birks' 
Law [25] to model recombination effects in LAr due to 
highly ionizing particles [26 . 

Simulated Monte Carlo (MC) single-monopole sam- 
ples are used to determine the efficiency as a function 
of the transverse kinetic energy E T ln = E kln sin 9 and 
pscudorapidity 77 for various monopole masses. For the 
Drell-Yan process, it is assumed that spin- 1/2 magnetic 
monopoles are produced in pairs from the initial pp state 
via quark-antiquark annihilation into a virtual photon. 
MadGraph [37] is used to model this process by assum- 
ing leading-order Drell-Yan heavy lepton pair production 
but making the replacement e — ¥ g/3 to reflect the mag- 
netic coupling in Eq. [2] In the Drell-Yan samples, the 
CTEQ6L1 [28 parton distribution functions are used and 
Pythia version 6.425 [39] is used for the hadronization 
and the underlying event. Only Drell-Yan monopoles 
with transverse momentum p^ > 200 GeV are processed 
by the simulation since lower monopoles fail to reach 
the calorimeter. For all the simulated samples, both the 
monopoles and the antimonopoles are assumed to be sta- 
ble and all final-state particles are processed by the sim- 
ulation of the ATLAS detector. Additional pp collisions 
in each event are simulated according to the distribution 
of pp interactions per bunch crossing in the selected data 
period. 

A simple algorithm is used to preselect events with 
monopole candidates for further study. Only events col- 
lected with a single-electron trigger with transverse en- 
ergy threshold Et > 60 GeV are considered. Monopoles 
that fulfill this energy requirement travel fast enough to 
satisfy the timing requirements of the trigger. Very high 
energy monopoles that exit the EM calorimeter are re- 
jected by a veto on hadronic energy that is intrinsic to 
the single-electron trigger. This trigger recorded an in- 
tegrated luminosity of 2.0 fb _1 during 2011 data-taking, 
defining the dataset used for this search. 

The reconstructed EM cluster is then required to have 
E T > 65 GeV and |?7| < 1.37. The trigger efficiency is 
independent of Et for Et > 65 GeV, motivating the for- 
mer requirement. The 77 requirement ensures that the 
EM cluster is in the barrel region of the LAr calorimeter, 
where the two-dimensional spatial resolution is uniform. 
If two or more EM clusters in an event satisfy these crite- 
ria, only the cluster with the highest energy is considered 
as a monopole candidate. 



In the barrel region, the monopole typically traverses 
35 TRT straws and its high ionization ensures that most 
of these register high-threshold hits. Furthermore, as 
each 5-ray produced by the monopole ionization deposits 
~ 2 keV in a straw, the combined energy deposited by 
multiple (5-rays crossing a single TRT straw gives rise to 
additional high-threshold hits. The large number of 6- 
rays bend in the 2 T magnetic field in the r-cj> plane; 
therefore, the monopole trajectory appears as a ~ 1-cm- 
wide swath of high-threshold TRT hits. The fraction of 
TRT hits that exceed the high threshold in the vicinity 
of the path of an ionizing particle is therefore a power- 
ful discriminator between the monopole signal and the 
background. The <f> position of the EM cluster is used to 
define a road of width A<fi = ±0.05 rad from the beam- 
line to the cluster. At least twenty high-threshold TRT 
hits must be present in the road. In addition, at least 
20% of the TRT hits in the road must be high-threshold 
hits. 

After the preselection, a more refined TRT hit count- 
ing algorithm is used to distinguish the signal from the 
backgrounds. A histogram with a bin width of 0.8 mrad 
is filled with the <j> distribution of the high-threshold hits 
in the previously defined road. The location of the high- 
est bin is used to calculate the center of a new road. In 
the TRT barrel, a rectangular road of ±4 mm in the r— <j> 
plane is used and the hits are counted. In the TRT end- 
cap, a wedge-shaped road of width A(j) = ±0.006 rad is 
used. These roads are wide enough to encompass two 
neighboring straws, taking into account the monopole 
trajectory and the associated (5-rays, but sufficiently nar- 
row to ensure that the fraction of hits that exceed the 
high threshold is insensitive to the presence of neighbor- 
ing tracks. In the barrel region, the number of hits in the 
road is required to be greater than 54. An //-dependent 
requirement on the number of hits in the road is applied 
in the endcap and barrel-endcap transition region to ac- 
count for the TRT geometry. 

Energy loss by bremsstrahlung and e + e~ pair produc- 
tion is negligible for magnetic monopoles in the mass 
and energy range considered herein. Therefore, mag- 
netic monopoles give rise to a narrow ionization energy 
deposit in the LAr calorimeter, the size of which pro- 
vides another powerful discriminator of the monopole 
signal from backgrounds such as electrons and photons, 
which induce an EM shower. The variable used is or, 
the energy-weighted two-dimensional rj-(f> cluster disper- 
sion in the second layer of the EM calorimeter, which 
has the highest two-dimensional spatial resolution. The 
dispersion <tr is calculated from the energies deposited 
in a 3 x 7 array of cells centered around the most en- 
ergetic cell of the EM cluster: ctr = + °f; 7 where 

g\ = E (EiStf) IY.E,, - [E (E^) /S^] 2 , 5^ is the de- 
viation in 4> between cell i and the most energetic cell 
and Ei is the energy of cell i; a'i is defined similarly. 
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FIG. 1. High-threshold TRT hit fraction, f HT , versus EM 
cluster dispersion, <tr. The circles represent 1000 simulated 
single monopoles with mass 800 GeV. The crosses represent 
ATLAS data. The regions marked A-D are discussed in the 
text. 



The fraction of TRT hits that exceed the high thresh- 
old in the narrow road pointing to the cluster, /ht, and 
the cluster dispersion, or, are thus chosen as the distin- 
guishing variables between the signal and background, 
and are shown in Fig. [l] The main physics background 
sources are high energy electrons, photons and jets, which 
exhibit no correlation between these variables in sim- 
ulated processes. The background and monopole MC 
samples are used to define an approximate signal re- 
gion. Toy MC events are generated by sampling the one- 
dimensional or, and /ht distributions for data outside 
this approximate signal region. The borders of the signal 
region are then tuned for optimal signal sensitivity by 
replacing the background MC events with these toy MC 
events. The final signal region A is defined by or < 0.017 
and /ht > 0.7. 

The efficiencies, which include trigger, reconstruction 
and selection effects, in the two-dimensional E^ ln versus 
7] plane are obtained from the simulated single-monopole 
samples. A fiducial region for each monopole mass is de- 
fined by the E^ ln range in which the efficiency is 0.80 
or higher in the \rf\ < 1.37 region. Figure [2] shows the 
efficiency versus -E T ln , averaged over \rj\ < 1.37. For 
monopoles with a mass of 200 GeV, the minimum trans- 
verse kinetic energy (E^ ln ) m i n where the efficiency rises 
above 0.80 is 700 GeV. For monopoles with a mass be- 
tween 500 GeV and 1500 GeV, (E$ m ) min is 600 GeV. 
Monopoles with lower E^ n fail to penetrate to the EM 
calorimeter and therefore do not satisfy the trigger re- 
quirements. Monopoles with very high E^ n exit the 
EM calorimeter and are rejected by the hadronic veto 
of the electron trigger. A common upper value of E^p n = 




E?" [GeV] 

FIG. 2. Efficiency versus -Br", averaged over \q\ < 1.37, for 
single monopoles of mass 200 GeV and mass 1500 GeV. 



1400 GeV is used for the fiducial region of all monopole 
masses. As the minimum efficiency is 0.80 in the fiducial 
region, a common value of 0.80 is used in the determina- 
tion of the upper cross-section limit. 

The efficiencies can be under- or over-estimated for sev- 
eral reasons. These effects are described below and the 
relative systematic uncertainties for each effect are given. 
1) Cross-talk in the second EM layer in the 4> direction is 
not modeled in the simulation. The energy is rcweightcd 
assuming 1.8% cross-talk [3U] and the cluster dispersion, 
or, recomputed. The efficiency is reduced and the re- 
sulting relative shift of —1.7% for single monopoles is 
taken as a one-sided uncertainty. 2) The simulation un- 
derestimates the TRT occupancy in the data by up to 
20%; therefore, the number of low-threshold hits (those 
unlikely to come from the monopole or related 5-rays) 
is increased by 20%. The resulting relative uncertainty 
is —1.3%. 3) The modification to Birks' Law is var- 
ied between its upper and lower systematic uncertain- 
ties [53], yielding a relative uncertainty of +1.8% and 
+1.5%, respectively. 4) The production of 5-rays is var- 
ied by 3% [T5] and the resulting uncertainty is negligible. 
5) The Geant4 "range cut" [23] controls the minimum 
kinetic energy threshold below which 5-rays are not prop- 
agated explicitly. This parameter is reduced from 50 /im 
to 25 /im in the TRT simulation. The resulting relative 
uncertainty is +0.14%. 6) The material in the inner de- 
tector, in the barrel cryostat and in between the cryostat 
and the first layer of the EM calorimeter is increased by 
5%, 10% and 5%, respectively, in the simulation [31]. The 
resulting —0.74% relative uncertainty is taken as sym- 
metric. Including an uncertainty of ~1.7% to account 
for the limited number of MC events, the total upper 
and lower relative uncertainties on the efficiency for sin- 
gle monopoles are +2.6% and —2.8%, respectively. 

The efficiencies for monopoles produced with Drell- 
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Yan kinematic distributions are given in Table [T] for each 
mass. Only masses up to 1200 GeV are considered, tak- 
ing into account the phase space limitations for pair pro- 
duction. The total relative uncertainties, which reflect 
the same systematic variations described above, are also 
given. The efficiencies and their associated systematic 
uncertainties reflect large losses due to acceptance, since 
many Drell-Yan monopoles have insufficient energy to 
reach the calorimeter. 



TABLE I. Efficiencies and their relative uncertainties in per- 
cent for Drell-Yan pair- produced monopoles of various masses. 



Mass (GeV) 


200 


500 


800 


1000 


1200 


Efficiency 


0.011 


0.048 


0.081 


0.095 


0.095 


Relative uncertainty 












Upper (%) 


+32 


+24 


+22 


+23 


+20 


Lower (%) 


-36 


-23 


-22 


-25 


-25 



The background in the signal region is predicted di- 
rectly from the data. The two-dimensional plane in Fig.[l] 
is divided into quadrants, one of which is dominated by 
signal (region A), and three others that are occupied 
mainly by background (regions B, C and D). The ra- 
tio of background events in signal region A to events in 
background region B is expected to be the same as the 
ratio of background events in regions C to D. This as- 
sumption is incorporated into a maximum likelihood fit 
to determine the estimated numbers of signal and back- 
ground events in signal region A. The inputs to the fit 
include the observed event yields in quadrants A through 
D, which are 0, 5, 16 and 7001, respectively, the efficien- 
cies and associated systematic uncertainties that have al- 
ready been discussed, and the integrated luminosity and 
its 3.7% uncertainty [32]. For each monopole mass, the 
rate of appearance of signal events in quadrants B and C, 
as predicted by the simulation, is also taken into account. 
According to the simulation, no signal event appears in 
quadrant D for any monopole masses. The fit predicts 
0.011 ± 0.007 background events in the signal region. 

Using the results of the maximum likelihood fit, the 
upper limits on the production cross sections at 95% con- 
fidence level are calculated using the profile likelihood ra- 
tio as a test statistic 33 . The results extracted using the 
CL S method [34] are shown in Fig. [3j Since the number 
of expected background events is very small and no event 
is observed in the signal region, only the observed limits 
are shown. The fluctuations of the observed limit in the 
fiducial region originate from variations of the nuisance 
parameters used in the profile likelihood ratio. 

In conclusion, the upper limit on the production cross 
section at 95% confidence level is 2 fb for Dirac mag- 
netic monopoles with the minimum unit magnetic charge 
and with mass between 200 GeV and 1500 GeV. No as- 
sumption about the production mechanism is made for 
this result, which is valid in the fiducial region given 



by \rj\ < 1.37 and (^ in ) mi „ < £ kin sin6l < 1400 GeV, 
where (-Erp m )roin is 600 GeV for monopoles with a mass 
between 500 GeV and 1500 GeV. For monopoles with 
a mass of 200 GeV, (£^ in ) min = 700 GeV. Assuming 
the kinematic distributions from Drell-Yan pair produc- 
tion of spin-1/2 Dirac magnetic monopoles, the upper 
limit on the cross section at 95% confidence level varies 
from 145 fb to 16 fb for monopoles with mass between 
200 GeV and 1200 GeV. These results extend the upper 
limits on the production cross-section for monopoles in 
this mass region established by preceding experiments. 
This collider search is the first that yields constraints on 
magnetic monopoles with masses greater than 900 GeV. 
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FIG. 3. Upper limits on the monopole production cross sec- 
tions at 95% confidence level. The solid line is the limit for 
single monopoles in the fiducial region and the dashed line is 
the limit assuming the kinematic distributions from Drell-Yan 
(DY) monopole pair production. 
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